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Abstract 

We establish new lower bounds for the convergence radius of the Mayer series and the 
Virial series of a continuous particle system interacting via a stable and tempered pair 
potential. Our bounds considerably improve those given by Penrose and Ruelle in 1963 for 
the Mayer series and by Lebowitz and Penrose in 1964 for the Virial series. To get our results 
we exploit the tree-graph identity given by Penrose in 1967 using a new partition scheme 
based on minimum spanning trees. 
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1 Introduction 

Continuous particle systems are an important subject of investigation in rigorous statistical mechan¬ 
ics. Since the thirties it was clear that the equation of state of a non ideal gas in thermodynamics 
can be deduced from statistical mechanics via the Mayer series of the pressure in the grand canon¬ 
ical ensemble. The Mayer series, proposed by Mayer (see, e.g., m and references therein) is an 
expansion in powers of the fugacity A of the logarithm of the Grand Canonical partition function 
of the the system under study, whose n*'^-order coefficient is given by a sum of terms indexed by 
the labeled connected graphs between n vertices. This series was at the time only formal and 
the question about its convergence uniformly in the volume where the gas is confined, remained 
an enigma during two decades. The difficulty in dealing with this problem was basically due to 
challenging combinatorial structure of the Mayer coefficients. 

In the sixties, Penrose m and independently Ruelle m proved the convergence of the Mayer 
series for a very large class of continuous particle systems. Namely those interacting via a stable 
and tempered pair potential. The bound on the convergence radius R of the Mayer series obtained 
in m and m still stands as the best available in the literature. It is worth to remind that the 
first singularity of Mayer series is in general not on the positive real fugacity axis (e.g. it is in the 
negative axis at A = — ii for repulsive potentials, see m) so that the value A = R for the fugacity is 
not directly related to any physical phase transition point. However a lower bound for R as sharp as 
possible is a physically relevant information since it allows to maximize the region of analyticity of 
the pressure of the system, i.e. where the system is surely in the gas phase and no phase transitions 
occur. In this respect, there have been only very few improvements on the Penrose-Ruelle bound 
for R by considering some restricted sub-classes of stable and tempered potentials. Basically such 
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improvements can be summarized as follows. Brydges and collaborators improved the Penrose- 
Ruelle bound for absolutely summable pair potential (see e.g. [H] and references therein) and 
Basuev gave in [2] an improvement for a significant class of potentials, introduced by himself in [T], 
which is sufficiently large to include Lennard-Jones type potentials (see also recently m)- 
The method used to get the general bounds of 1963 were based on the so-called Kirkwood-Salzburg 
equations (see [TS] and references therein). It is an indirect method which aims to control simulta¬ 
neously all correlation functions of the systems by analyzing an infinite set of equations involving 
them. On the other hand, the methods used to get improvements cited above for subclasses of sta¬ 
ble and tempered potentials are all based on trying to bound directly the Mayer coefficients. The 
main technical tools to obtain this task are the so-called tree-graph identities. Several alternative 
identities have been proposed along the last four decades (see mm and references therein). 

The very first of these tree-graph identities, proposed by Penrose m in 1967, was based on the 
existence of a map M (a so-called partition scheme) from the set Tn of the labeled trees with n 
vertices to the set Qn of the labeled connected graphs with n vertices. This map is able to induce a 
partition of the set Qn into blocks indexed by the elements of set Tn (details in Section [3|). Penrose 
gave an explicit construction of such a map and used it to rewrite the Mayer coefficients as sum 
of terms indexed by trees rather than by connected graphs. Using this alternative expression of 
the Mayer coefficients, Penrose was then able to bound them directly, reobtaining the same bound 
of 1963 but at the cost to impose a further condition on the pair potential, beyond stability and 
temperedness. Namely, the pair potential needed to have an hard-core. This fact, we think, has 
led researchers to believe that the identity proposed by Penrose was useful only for systems with 
hard-core and in fact in the literature it is mentioned only for such kind of systems (see, e.g., [3 [20] 
and reference therein). 

In the present paper we disclaim this belief and show that the Penrose tree-graph identity can be 
used to strongly improve the old Penrose-Ruelle bounds for general stable and tempered potentials. 
The main idea is based on the definition of a new map M, different from the one originally proposed 
by Penrose, which still is able to make a partition of the set of connected graphs whose blocks are 
indexed by minimum spanning trees, where the minimality is basically on the energies of the edges 
of the tree. Such new partition scheme, as we will see, allows to use in an optimal way the stability 
condition of the pair potential. 

The rest of the paper is organized as follows. In Section [5] we remind the basic concepts for 
continuous systems of classical particles, we give the definitions of stability and temperedness and 
write down the Mayer series explicitly. We end the section by recalling the Penrose-Ruelle result 
(Theorem 0) and by stating the main results of the paper. Namely, in Theorem [1] we present 
the new bound on the n*^-order Mayer coefficient and consequently the new lower bound on the 
convergence radius of the Mayer series and in Proposition [1] we give the main technical tool to proof 
Theorem (T] id est a new tree graph inequality, deduced from the Penrose tree-graph identity. In 
Section 3 we prove Proposition [TJ We start by recalling the Penrose identity for a generic partition 
scheme M, we then introduce our new partition scheme and show how to use it to obtain the 
inequality given in Proposition [TJ In Section 4 we conclude the proof of Theorem [T] Finally, in 
Section 5 we give a comparison with previous results and add some concluding remarks. 
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2 Model and Results 


2.1 Model 

We consider a continuous system of classical particles in the d-dimensional Euclidean space 
Denote by Xi G the position vector of the particle of the system and by HxiH its Euclidean 
norm. We suppose hereafter that particles interact through a pair potential V : —>■ M U {-|-oo} 

such that V(—x) = V(x) and we set shortly R* = M U {-|-oo}. Given a configuration (xi,..., Xn) G 
^dn q£ system such that n particles are present, the energy U (xi,..., Xn) of this configuration 
is dehned as 

U{xi,. . . ,Xn) = ^ V{Xi-Xj) 

We work in the Grand Ganonical Ensemble, where the statistical mechanics of the system is gov¬ 
erned by the following partition function 


OO , ^ « « 

Ea(A,/3) = V — / dxi... / 


( 2 . 1 ) 


where A C R'^ is typically a d-dimensional cube with center in the origin, A G (0, -|-oo) is the activity 
of the system and j3 G (0, -|-oo) is the inverse temperature. 

The connection with thermodynamics is obtained by taking the logarithm of the partition func¬ 
tion. Namely, the pressure P\{(3,X) and the density /9 a(/ 3, A) of the system at fixed values of the 
thermodynamic parameters inverse temperature /3, fugacity A and volume A, are given by 

^A(/3,A) = ^logHA(A,/3) (2.2) 

dA(/3,A) = A|.iogHA(A,/3) (2.3) 

It is a long known fact [12] that logHA(A,/3) can be expanded as a power series of the fugacity. 

Namely, 

.. CX> 

— logHA(A,/l) = X + J2Cn{/3,A)X^ (2.4) 

The series above is the so-called Mayer series and its coefficients Cn(/3,A), known nowadays as 
Mayer coefficients, admit the following explicit expression. 


CnW,A) = 


|A| ra! A 


dxi 


dXr 


E n 

g&Qn {i,j}€Eg 


^-pV{xi-Xj) 


- 1 


(2.5) 


where Qn denotes the set of all connected graphs with vertex set {l,2,...,n} (a.k.a. labeled 
connected graphs with n vertices) and Eg denotes the edge set of <7 G 

Clearly, via Mayer series one can write both pressure (12.2p and density (12.31) as expansions in power 
of the fugacity A. 

Some conditions on the potential V must be imposed to have hope to control the series. Stability 
and temperedness are commonly considered as minimal conditions to guarantee a good statistical 
mechanics behavior of the system (see, e.g., [T8] and El)- 
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Definition 1 (Stability) A pair potential V is said to he stable if 


B := sup sup - U{xi,... ,Xn) <+00 

»^>2 (a;i,...,Xn)GR‘'" 

We call B the stability constant of the pair potential V. 


( 2 . 6 ) 


Definition 2 (Temperedness) A pair potential V is said to be tempered if there exists ro > 0 
such that 


||x||>ro 


\V{x) \ dx < oo 


(2.7) 


Observe that the stability constant of any tempered potential V is non-negative. 

As mentioned in the introduction, the best rigorous upper bound on |C'n(/3, A)| so far (and hence 
the best lower bound on the convergence radius of the Mayer series) for stable and tempered pair 
potentials is that obtained by Penrose and Ruelle in 1963 [151119] . 

Theorem 0 (Penrose Ruelle 1963) Let V be a stable and tempered pair potential with stability 
constant B. Then the n-order Mayer coefficient Cn{/3,A) defined in 112.5]) is bounded by 


n\ 


n—l 


( 2 . 8 ) 


where 


C{fi) = [ dx - 1| 


(2.9) 


Therefore the Mayer series (2-4) converges absolutely, uniformly in A, for any complex A inside 
the disk 


|A|< 


1 


e2/3s+iC(/3) 

I.e. the convergence radius R of the Mayer series admits the following lower bound 


( 2 . 10 ) 


R > RpK 


1 

e‘^dB+ic{i3) 


( 2 . 11 ) 


Observe that the non-negative quantity C{I3) defined in (j2.9j) is finite if V is stable and tempered 
(see e.g. [I8]). 

2.2 Results 

As said in the introduction, the bound on the absolute value of the n**^-order Mayer coefficient 
C'n(/3, A) appearing in Theorem 0 has been obtained “indirectly” using the method of Kirkwood- 
Salzburg equations. Our main result is the following Theorem [1] below which improves strongly 
Theorem 0 via a direct bound on Cn{fi,A) by using the old Penrose tree-graph identity proposed 

in [T6] . 


4 








Theorem 1 Let V he a stable and tempered pair potential with stability constant B. Then the 
n-order Mayer coefficient Cn{/3ffi) is bounded by 


|C„(/3,A)| < 

n\ 


72—1 


where 


CiP) = [ 


1 _ e-/3|^(*)l 


dx 


( 2 . 12 ) 


(2.13) 


Therefore the Mayer series converges absolutely, uniformly in A, for any complex A inside the disk 


|A|< 


1 


ehB+ic{ffi 

I.e. the eonvergence radius R of the Mayer series admits the following lower bound 

R> R* = ^ 


(2.14) 


ehB+ic{ffi 


(2.15) 


Observe that, according to the definitions (j2.9p and (I2.13p . we have 0 < C{f3) < C{f3). 

Remark 1. The improvement given by Theorem [1] respect to Theorem 0 is twofold. First, the 
factor in (j2.15l) replaces the factor in (|2.1ip . Second, the factor C{[d) in (I2.15P replaces 

the factor C{/3) in (j2.1ip and clearly, recalling their definitions (I2.13P and (12.9p respectively, we 
have C'(/3) < C{(5) where the equality only holds if V is non-negative (purely repulsive). Moreover 
observe that while C{(5) grows at most linearly in j3, the factor C{I3) can grow exponentially with /3 
(for V with an attractive part). So, the ratio R*/Rp^i always greater than one, is the product of two 
factors, e^^ and [C'(/3)/C'(/3)], both growing exponentially fast with /3 when V has a negative part. 
Due to this exponential behavior in /3, our improvement on the Penrose-Ruelle bound is rather 
astonishing for attractive potentials (i.e. those with strictly positive stability constant B) at non 
small /3 (i.e. low temperatures). To give an idea, for the Lennard-Jones gas at inverse temperature 
(3 = 1, using the value Bpj = 8.61 for its stability constant (see M), our lower bound is at least 
8.5 X 10^ larger than the Penrose-Ruelle lower bound, while for /3 = 10 is at least 7.26 x 10^^ larger 
than the Penrose-Ruelle lower bound. On the other hand, for (3 small the exponential e^^ is near 
to one and C{(3) becomes very close to C((3) so that our improvement is much less sensitive at high 
temperatures. Finally, for purely repulsive pair potentials (i.e. those with stability constant B = 0) 
our bound coincides with the Penrose-Ruelle bound since in this case = 1 and C{(3) = C{(3). 

Remark 2. As mentioned above, the density of the system pj\^{(3,\) given in (12.31) can also be 
written, via (12.4p . in terms of a power series of the fugacity A. This series allows to express the 
fugacity A as a power series of the density p\ which, plugged into (12.41) . gives the so-called Virial 
series (see, e.g., Ha), i.e., the pressure as a function of the density and temperature, or, in other 
words, the equation of state of the system. In 1964 Lebowitz and Penrose [8] (provided the following 
lower bound on the convergence radius TZ of the Virial series of a gas of particles interacting via a 
stable and tempered pair potential. 


TZ > T^lp 


(apBc{j3) 


(2.16) 
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with 


[(1 + u)e ^ — l]t(; 


9{y) 


max 
0<-u)<l u 


The authors deduced such a bound, via Lagrange inversion, from the Penrose-Ruelle bounds (12.81) of 
the Mayer coefficients (see also [7] for an alternative method that deal directly with the coefficients of 
the Viral series written in terms of two-connected graphs). If one redoes the calculations performed 
by Lebowitz and Penrose using the new upper bound (12.1211 for the absolute value of the n^'^-order 
Mayer coefficient (with C{/3) in place of C{/3) and in place of then one immediately 

obtains the following lower bound for the convergence radius TZ of the Viral series 


7^ > 


g(e^^) 

ef^^C{l3) 


(2.17) 


Observing that 0.1448 < g{u) < e~^ for u G [l,oo) (see, e.g.. Remark 2 below Theorem 1 in |13jL 
we have that the ratio 77.*/77 lp between the new lower bound (I2.17P for the convergence radius 
of the Viral series and the old Lebowitz-Penrose bound given in (|2.16|) is of the same order of 
magnitude of the ratio R*/Rpr calculated above. 

The proof of Theorem [1] relies on a non trivial tree-graph inequality on the absolute value of the 
integrand of the r.h.s. of (|2.5p (the so-called Ursell functions) as far as stable pair potentials are 
considered. This new tree-graph inequality, an additional result of the present paper which we 
consider interesting per se, can be stated as follows. 


Proposition 1 Let V be a stable pair potential with stability constant B. Then, for any n > 2 and 
any (xi,..., Xn) G the following inequality holds 


^ n - 1)1 < (2.18) 

g&Gn r&Tn {i,j}eE-r 


where Tn denotes the set of all trees with vertex set [n] (a.k.a. labeled trees with n vertices). 
The proofs of Proposition [1] and Theorem [1] are given in sections 3 and 4 respectively. 


3 Proof of Proposition [T] 

In this section we set [n] := {1, ...,n} and '■= {{i,j} ■ i,j G [n], i < j}. We recall also that Qn 
is the set of all connected graphs with vertex set [n] (i.e labeled connected graphs with n vertices) 
and Tn is the set of all trees with vertex set [n] (i.e labeled trees with n vertices). Given g G Gn, we 
denote Eg its edge set. Given g € Gn and g' G Gn such that Eg C Eg/, we say that (7 is a subgraph 
of g' and write g C g'. Note that the pair {Gn-, C) is a partially ordered set. Given g,g' G Gn such 
that g C g', we denote [g,g'] = {g" ^ Gn ■ 9 ^ 9 " G g'y. In other words [g,g'] is an interval in 
{Gn, C). Given a set S, we denote by |5| its cardinality. 
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3.1 General Penrose identity 

Penrose identity is a rewriting of the Ursell functions 




E n 

g&Gn {i,j}GEg 


Q-PV{xi-Xj) _ 2 


based on the existence of a map {partition scheme) from the set Tn of the trees with vertex set [n] 
to the set Qn of the connected graphs with vertex set [n]. 

Definition 3 A map M '■ Tn ^ Gn is called a partition scheme in the set of the connected graphs 
Gn if, for all T G Tn, T C M{t) (i.e. r is subgraph of M{t)) and Gn = Utg?;, [e-^('^)] where IjJ 
means disjoint union and [t,M{t)] = {g G Gn ■ x C g C M{t)} is an interval in Gn (with respect 
to the set-inclusion). 

Once a partition scheme in Gn has been given, we have the following identity 


Theorem 2 (General Penrose identity) Let V{x) be a pair potential. Let n > 2. Let M : 
Tn ^ Gn be a partition scheme in Gn- Then, for any {xi, ..., x„) G the following identity holds 


geGn {i,j}GEg 


^-pVixi-Xj) _ 

= E 



r&Tn 



-0 E V(xi-Xj) 


(3.1) 


Proof . Let us pose shortly Vij = (dV{xi — xj). Since Gn is the disjoint union Gn = 1±Itg7;i ^ 
we can write 


E n 

g&Gn {i,j}&Eg 


- 11 = 


1) = E E n 1) = 

r€Tn gG[T,M(r)] {i,j}eEg 


E n 

reTn {i,j}€ET 

E n 

rGTii {i,j}€ET 


(e-^- - 


fe-^- - 


1 ) E n = 

E CE^CE Gd}<^Eg\Er 

1 ) E n E""-!) = 

ECi?A<r(T)\®T {i,j}&E 


E 

tSTu 


n ( 

{i,j}&E. 


-Vii 

e 


-1) 


n - 1 )+ 1 ] 

{*d}G®M(T)\®T 


E 

r&Tn 


n (e - 1) 

{*d}GEr 


n 

{*d}€EAf(T)\ET 


-Vii 

e 


□ 


From Theorem [2] the corollary below easily follows. 
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Corollary 1 Let V{x) he a pair potential. Let n > 2. Let M ■ Tn ^ Gn be a partition scheme in 
Gn- Then, for any (xi,... ,Xn) G the following inequality holds 


g£Gn {i,j}£Eg 


reTn 


-h E V(xi-Xj) 


— e 


-P\V{xi-Xj)\ 


where Eif denotes the set of edges of the tree r with non-negative energy. That is, 

= {{i,j} e Er : V{xi- Xj) > 0} 

Proof. Let us pose once again Vij = f5V{xi — Xj). By Theorem [2] we have that 


E n 

geSn {i,j}&Eg 


T&Tn {i,j}&ET 


Vii 


Using now the trick proposed in m we observe that, for any t ^Tn 


n - ii=[ n 

{i,j}&Er {i,j}&Er 


— e 


-Toh 




(3.2) 

(3.3) 


(3.4) 


so that 


|e - l| 



_{i,j}SEr- 


{i,j}eEr 




(3.5) 


Inserting now (13.5p into ()3.4p we get (13.2p . □ 


3.2 Partition scheme via minimum spanning tree 

The key point now is to find a partition scheme M in such a way that is possible to find a good 
bound for the factor exp{—/3 Yl{i j}£EM( )\E''' “ ®i)} appearing in the r.h.s. of (13.2p . 

We will construct explicitly our partition scheme M ■ Tn ^ Gn first defining an auxiliary map 
T ■ Gn ^ Tn and then deriving M from T according to the following proposition. 


Proposition 2 The following statements are equivalent. 
1. There are two maps 


T 



such that T ^(r) = {g € Gn ■ t C g C M{t)} for every r G 7^. 


2. M is a partition scheme in Gn- 
















Proof. 1 => 2. Since g G T~^{T{g)), we have T{g) C g for all g G Qn- In particular, for every 
tree r we have T(t) C r which implies T(t) = r because both are trees. I.e., T is surjective and 
thus the intervals T~^(t) are nonempty. This implies that Qn is the disjoint union of the intervals 

M{t)]. Hence in view of Dehnition [3] we conclude that M is a 

partition scheme in Q^. 

2 1. If M is a partition scheme then for any g & Qn there exists a unique tree t a Tn such that 

g G Therefore we can define the map T from Qn to Tn such that, for all g G 

T{g) = T. □ 

We thus start by first defining the map T from Qn to Tn- In order to do that, assume we have 
a function defined on the edges of the complete graph, / : En M*. Then, for every connected 
graph g & Qn there is at least a tree t C g, among the trees t' C g, which minimizes the value 
of / /(®)- This tree is called a minimum spanning tree of g w.r.t. /. If / is such that this 

minimum spanning tree is unique for each graph g (z Qn, then / induces a map, say Tf, from Qn to 
Tn which associates to each g Qn this unique minimum spanning tree. 

Let now H be a stable and tempered pair potential in and let (xi,..., x„) G be given. We 
would like to consider as the function / : En M* the one that associates to each {i,j} G En the 
value V{xi — xj). The problem is that such /, depending on the potential V and the conhguration 
{xi, ... ,Xn) G chosen, in general does not guarantee that each g ^ Qn contains a unique 
minimum spanning tree w.r.t. to /. In order to avoid multiple minima, we will therefore modify 
the codomain of the function / by adding n(n— l )/2 auxiliary “coordinates” (which will take only 
integers values), each of them one-to-one associated to an edge of En- These coordinates will permit 
to distinguish eventual multiple minimum trees. 

To carry out this procedure, we need to widen our framework and consider, instead of functions 
from En to M*, a suitable class of more general functions / : En —?> IK with IK being a totally ordered 
Abelian monoid and having the key property that YIegEt /(®) different for different trees t € Tn- 
We recall that a totally ordered Abelian monoid (shortly tomonid) is a structure (IK,-|-,0, >) such 
that (IK,-l-,0) is an Abelian (i.e. commutative) monoid, (IK, >) is a totally ordered set and for 
all x,y,z G IK we have that x > y implies x + z > y + z (i.e. the total order > is translational 
invariant). 

Definition 4 Let f : En —>■ IK where IK zs a totally ordered Abelian monoid. We say that f is 
admissible if for any t,t' GTn such that r 7 ^ r' rce have that Ye€E{T) /(®) Yle&Ep') /(®)- 

Once an admissible function / : En IK has been given, we can define the maps Tf from Qn to Tn 
Mf from Tn to Qn as follows. 

Definition 5 Let IK be a totally ordered Abelian monoid. Let f : En be an admissible function. 
Then, for every g G Qn there is a unique spanning tree t C g for which Ye&Ep) /(^) minimum. 
We define the map Tf : Qn ^ Tn such that Tf{g) is this unique minimum spanning tree of g. 

Definition 6 Let IK 6 e a totally ordered Abelian monoid and let f : En -G- ^ be an admissible 
function. We define the map Mf : Tn ^ Qn such that Mf^r) is the graph on the vertices [n] whose 
edges are the {z, j} such that f{{i,j}) > fie) for every edge e G E^ belonging to the path from i to 
j through r. 
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Thus we have constructed 


Tf 

Qn < " Tn 

Mf 

Observe that r C Mf{T) and Tf{g) C g. Moreover, these two maps Tf, Mf satisfy the hypothesis 
of the Proposition [2j This is the content of the following lemma. 

Lemma 1 Let f : En —>■ K 6e an admissible function and r G Tn- Let Tf and Mf be the maps 
given in Definitions\^and\^ respectively. Then 

Tf^ir) = {g egn- T C g C 

and therefore Mf is a partition scheme in Qn- 

Proof. Let g G We have r = Tf{g) C g. Now take {i,j} G Eg, and let e G Sr be any 

edge belonging to the path from i to j in r. Consider t' the graph obtained from r after replacing 
the edge e by {i,j}. Clearly t' is connected and has n — 1 edges, so it is a tree. By minimality of 
r we must have f{d) < f{{i,j}), whence {i,j} G S^^(r)- Therefore g C Mf{T). 

Conversely, let r C 5 C Mf^r). We must show Tf{g) = r. By cardinality, it suffices to show 
Tf{g) C r. Proceeding by contradiction, take {i,j} G \ E^- Consider the path p'^{{i,j}) in 

r joining i with j. Since Tf{g) C Mf(T), f({i,j}) is greater or equal than the corresponding value 
for any edge in the path p'^({i,j}). If we remove {i,j} from Tf{g), the tree splits into two trees. 
Necessarily, at least one of the edges in the path p^{{i,j}) joins a vertex of one tree with a vertex 
of the other. Thus, by adding this edge we obtain a connected graph with n — 1 edges, a new tree, 
which contradicts the minimality of Tf{g). 

We now provide an explicit construction for the admissible function / of Definition [5l 
We start by specifying the totally ordered Abelian monoid K. Assume to have chosen an order 
in En (e.g. the lexicographic order: {i, j} < {i'ij'f if either i < i! or i = i! and j < /). Let 
No = N U {0} and consider the set 


\En\ times 

N®" = No X X No 

such that the entry (with 1 < m < \En\) of an element x G N^"' corresponds to the edge 
w.r.t. the order chosen of En- Now we set 

K = M* X (3.6) 

An element of K is thus an ordered {\En\ + l)-tuple such that the first entry is a real number while 
the remaining \En\ entries are natural numbers or zero. 

The set K defined in (13.61) has a canonical structure of Abelian monoid and can also be endowed 
with a natural total order by considering the lexicographical order on M* x N™, prioritizing the 
entries from left to right. It is easy to show that such an order is translational invariant respect to 
the standard sum in K which is thus a totally ordered Abelian monoid. 

We are now ready to define the function / : En —>■ IK. Let l{ij} be the element of N^" with all 
entries zero except the one at the position corresponding to the edge {i,j} which is equal to one. 
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Definition 7 Given a pair potential V and given (xi,..., x„) G define the function 

f: En^K: {i, j} - Xj) x l{jj} (3.7) 

Remark. The function / above is, for every pair potential V and for any (xi,...,x„) G 
admissible according to Definition [5] and thus, according to Definition [6] and Lemma [U the map 
Mf : Tn ^ Gn a. partition scheme in Qn- This is onr new partition scheme. The advantage 
to consider this new partition Mf scheme instead of the one originally proposed by Penrose is 
manifestly clear in the following key lemma which practically concludes the proof of Proposition [TJ 

Lemma 2 Let P : — >■ M* 6e a stable pair potential with stability constant B, and t G Tn- Then, 

for every {xi,...,Xn) G 


M f{T)\Et 


Xj) > —Bn 


(3.8) 


where f is the function given in Definitions^ Mf{T) is the graph given in Definition\^ and Efi is 
the subset of Er defined in liS. 3\) . 

Proof. In the following, Vij denotes V{xi — xj), and we will make implicit use of the following 
trivial fact: 

(x, a) > (x', fj') X > x' for (x, a), {x , a') G M* x 

Now we proceed to show that the inequality ()3.8p holds true. The set of edges E-j- \ Efi forms the 
forest {ri,...,rfc}. Let us denote Vr„ the vertex set of the tree Tg of the forest. Assume i G Vr^, 
j G Vrf If a 7 ^ 6, the path from i to j through r involves an edge e in Efi. Thus, if in addition 
{Lj} G L'M/(r )5 have Vij > Pe > 0. If a = 6, the path from i to j through r is contained in Ta- 
Thus, if in addition {i,j} ^ we must have pj < Pe < 0 for some edge e in that path. This 

allows to bound: 


k k 

Y ^L>Y Y yL-^Y-^^rfiBy-nB □ 

^=1 PJlCV., S=1 


3.3 Conclusion of the proof of Proposition [T] 

We now have all ingredients to conclude the proof of Proposition [TJ Indeed, for any fixed configu¬ 
ration (xi ,... ,Xn) G by Corollary [1] and Lemma [1] we have 


Y n - 1) 

< 

(1 _ Q-h\V(xi-Xj)\^^ 

g&Gn {i,j}^Eg 




-hE 




+ V(xi-Xj) 


< 


E n 

T&Tn {iJ}GT 

where / is the admissible function given in (|3.7I) . Mj is the partition scheme given in Definition [6| 
and where to get the inequality of the last line we have used Lemma [2j 
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4 Proof of Theorem [T] 


To deduce Theorem [T] from Proposition [T] we just need to show the following easy lemma. 
Lemma 3 For any t ^Tn it holds 


dxi ... dx, 


n 

{i 


1 - < |A| [(7(/3) 


n—1 


Proof. Without loss of generality we can always assume (eventually through a renomination of the 
indices in the integral of the l.h.s. of equation above) that E^- = {2, j 2 }, { 3 ,^ 3 },... {n,jn} with 
jk < A; for all A: = 2,..., n, so that 




n 






dXn 

k=2 


_ Q-h\V{x,,-Xji^)\ 


Define the following change of variables (xi,..., Xn) —>■ (yi,..., Vn) in the integral 

yi=xi, yk = Xk-Xj^, \/k = 2,...,n (4.1) 

The Jacobian matrix of such transformation (j4.ip . being lower triangular with entries of the diagonal 
all equal to one, has determinant equal to one. Therefore we have 


< 


[ dxi... [ dXn Tf < 

Jk Jk ^ ^ 

dyi dy 2 ... dynT] 

Jk Jr3 Jk3 V 




= |A| 


(X - 


1=2 
-| n—1 


= |A| C{/3) 


n—1 


□ 


Theorem [T] follows now straightforwardly. Indeed, from Proposition [T] and Lemma [3] we get the 
following bound for the absolute value of the Mayer coefficient (^^(/S, A). 


< 


1 1 
|A| n! 

g/3Bn X 

|A| n 

^l3Bn X 


/ dxi . 

. / dXn 

E 

n [ 

Ik 

Ja 

g^Gn 

{hj}&Eg 


- f dxi ■■■ I dXn E If (1 

J TV J a. f • 'I y— rp 


< 


_ Q-d\V{xi-Xj)\^ _ 


|A| 


r&Tn {i,j}€ET 
[ dxi... [ dxn n (1 - < 




tGTu 


n! L 


= 


n 


n—2 


ni 


c{P) 


1 n—1 


where in the last line we have also used the Cayley formula 1 = |7^| = n” ^ (see [1]). This 

concludes the proof of Theorem [TJ 
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5 Comparison with recent results and concluding remarks 


As observed in Section 2, the improvement of Theorem[T]respect to the old Penrose-Ruelle bound is 
manifestly evident, as far as general stable and tempered pair potentials are concerned. There have 
been however recent works [laiiniin] which also obtain improvements on the convergence radius 
of the Mayer series for systems of particles interacting via a restricted class of stable and tempered 
pair potentials. In particular, in a recent paper m, which is a development of an early work by 
Basuev [5] , the authors obtain new bounds for the Mayer series convergence radius which are better 
than any precedent bound given in the literature (see in m the new bounds of Theorem 5, formula 
(4.5)), as far as the particles in the system interact via a so-called Basuev potential (see Definition 
3 in m)- This is a class of potentials which includes the classical case of the Lennard-Jones type. 
It is worth to say that, at least for the particular case of the Lennard-Jones potential, the optimal 
bounds presented in [TI] (obtained also using a recent result by one of us m) are slightly better 
than those obtained here. 

Indeed, as far as the Lennard-Jones potential V{r) = is concerned, the calculations given 

in Sec 5.3 of [n] show that, at /3 = 1 the lower bound for the convergence radius of the Mayer 
series of the Lennard-Jones gas is surely greater than 

7A 

where Blj is the Lennard Jones stability constant. The same calculation using instead the bounds 
of Theorem [1] yields for the same convergence radius a lower bound surely smaller than 

e“(^Lj+^) 

8J)8 

The results of this paper, as well as those given in laimiiiaiii], show that direct methods based on 
tree-graph identities are, in the end, much more effective than the old indirect methods of [laiii] 
based on K-S equations. On the other hand, the fact that the bound obtained in m beats our 
bound of Theorem [1] for some particular cases of stable and tempered pair potentials raises the 
question of whether it is possible to further improve the bounds of Theorem [1] in the general case. 
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